Rotating Spin-1 Bose Clusters 
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We propose a simple scheme for generating rotating atomic clusters in an optical lattice which 
produces states with quantum Hall and spin liquid properties. As the rotation frequencies increases, 
the ground state of a rotating cluster of spin-1 Bose atoms undergoes a sequence of (spin and orbit) 
transitions, which terminates at an angular momentum L* substantially lower than that of the boson 
Laughlin state. The spin-orbit correlations reflect 'fermionization' of bosons facilitated by their spin 
degrees of freedom. We also show that the density of an expanding group of clusters has a scaling 
form which reveals the quantum Hall and the spin structure of a single cluster. 



In this paper, we would like to point out some remark- 
able features of rotating clusters of spin-1 Bose atoms 
which result from the coupling of quantum Hall physics 
to spin symmetry through Bose statistics. Our studies 
are prompted by the recent experiments on rotating Bose 
gases and Bose Hubbard systems ffl , as well as the the- 
oretical realization of a 'fragmented' condensate in spin-1 
Bose gas J3|,|| , and the theoretical demonstration that a 
three dimensional Bose system can approach the quan- 
tum Hall limit at sufficiently large angular momenta || . 
As of now, one can release the spin degrees of freedom of 
a Bose gas in an optical trap ||, deposit a large amount 
of angular momentum to a Bose condensate m , and di- 
vide a Bose condensate into thousands of isolated clus- 
ters each with a few particles It is natural to ask 
whether these capabilities can be combined to produce 
a large collection of rotating and spin carrying clusters. 
In fact, such possibilities are being considered in some 
laboratories [?]. 

There are many reasons for studying these Bose clus- 
ters: Firstly, they are the systems where two types of 
strongly correlated ground states can be realized: quan- 
tum Hall droplets and fragmented Bose spin liquids. 
While the former is familiar to condensed matter physi- 
cists, the latter is new, with properties entirely different 
from the familiar single condensate (or coherent) states 
||. It is well known, however, that these states are very 
fragile as number of particle increases, easily giving way 
to the highly robust coherent states On the other 

hand, such fragility disappears in small clusters, and 
the system becomes strongly correlated in this regime. 
Secondly, quantum clusters are amenable to exact treat- 
ments, which may reveal physical principles operative in 
larger scale. Thirdly, amid the growing interest in using 
bosons in lattices to process quantum information, the 
study of quantum clusters remains a basic step. 

The results we will present are: a simple method to 
generate rotating clusters, the general features of their 
ground states, and the scaling behavior of the density 
of an expanding group of clusters, whose specific form 
reflects spin liquid and quantum Hall properties. Before 
proceeding, we first summarize the general features: 



(A): As the angular momentum L of the system in- 
creases, the interaction energy will be 'quenched' to zero 
at a critical value L* considerably lower than that of the 
boson Laughlin state. (B): As L increases from to L*, 
the total spin of the ground state is specified by its or- 
bital angular momentum, even though the interaction is 
invariant under separate spin and orbital rotation. (C): 
Both (A) and (B) reflect the phenomenon of 'fermion- 
ization' facilitated by the internal degrees of the freedom 
(i.e. spin) of the system. This is achieved by populating 
bosons in different spin-orbit states, a process automatic 
for fermions. (D): The energetics leading to (A) and (B) 
also implies that as the rotational frequency f2 increases 
toward the trap frequency w, the angular momentum of 
the ground state will increase in discrete steps until it 
reaches L*, then remains constant until the system be- 
comes unstable at ft > u>. We shall now derive these 
results. 

I. Rotating quantum clusters: their genera- 
tion and Hamiltonian: Consider adding a rotating 
quadrupolar potential W(r;t) — \{r\{t) — rf(i)), where 
r,(t) = v ■ Gi{t), ei(t) = xcosfif + ysinfM, e 2 (i) = 
ycosfH — xsinf^, and A is a constant. Such a poten- 
tial can be generated by a pair of off-center rotating 
laser beam as in the Paris experiment Q. The single 
particle Hamiltonian is then h(t) = T + U + W(r;t), 
where T = p 2 /2M, and U is the periodic potential of 
the optical lattice. Deep in the Mott limit, U is an ar- 
ray of deep wells at lattice sites {R}; the bottom of each 
is harmonic, U(r) — > l\ ar (x) = |Mw 2 x 2 as r — > R, 
r = x + R. These wells are so deep that atoms in dif- 
ferent wells are isolated from each other . Shifting the 
origin to R, the single particle Schrodinger equation is 
ihdt\^t) = [f+ L>har(x) + W(x + R;i)]|tf t ). Defining 
\%)= e -i[a(i)-P+b(i)-x]/n e -ifl(t) e -intLz | $( ^ and with 

appropriate choices of a, b and 9 p0[ , the Schrodinger 
equation can be written as ihdt\$>t) — fC\^t), where 
JC = T+Uhar(x-)+W(x; 0) — ilL z is time independent, and 
L z = xp y — yp x - Since W(x, 0) breaks cylindrical sym- 
metry, it will impart angular momentum to the system. 
Once the system acquires angular momentum, cylindri- 
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cal symmetry can be restored by turning off W, and K, 
reduces to the simple form JC = p 2 /2M + \Mu?y? -QL Z . 

Including particle interactions, the Hamiltonian of a 
rotating spin-1 Bose cluster is, in first quantized form, 
H - QL Z = t + V, where t = ^ V = E-*-^ # 



V; 



(c + c 2 F 4 • F,-) 6 (r, 



r^), where c and c 2 are 



density-density and spin-spin interactions [0,12j. The- 
oretical estimates show that c > for both ^Na and 
87 Rb; and that c 2 > (antiferromagnetic) for 23 Na and 
c 2 < (ferromagnetic) for 87 Rb (ll). Because of the 
similarity of scattering lengths in different angular mo- 
mentum channels, it is found that |c 2 |/c ~ 5% for both 
cases. From now on, we shall focus on values of Co and c 2 
appropriate for 23 Na and 87 Rb. For simplicity, we first 
consider zero magnetic field. Magnetic field effects will 
be considered at the end. 

The eigenvalues of K. are || E rhm . nz /h = (u + 
Q)n + (u> — f2)m + wn z , where n, m, and n z are non- 
negative integers. For very tight traps, Hui > Nc /a 3 
with a = yjTi/Muj, only states in the 'lowest Landau 
level' (LLL), (n = 0, m; n z — 0), will appear in the 
ground states. The wavefunction of the state (0, m; 0) 
is w m (x,y)f (z), where w m (x,y) = u m e~\ u 



fo(z) = e-^/^/iira 2 ) 1 /*, 
many body state is then 



I" 1 1 'v / 7ra 2 m!, 
and u — (x + iy)/a. The 
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<&> = B([u],[fA) II^( r 0|0> 
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where /(-) = / Ilili [ dr i /o(^) 



,-\ui\ 2 /2 
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(u\, U2, ■■■un), fM = 1,0,-1 labels the spin of the i- 
th boson, = (fix, fj,2, .../xat), -B([m],[/x]) is a sym- 
metric homogeneous polynomial of [u] labelled by [//], 
ip( r )fj, is the field operator in the LLL, ^( r ) = 
E m Wm(a;, y)fo{z)a\ ntl , and aJ„ M is the creation operator 
for the state (0, m; 0). In LLL, we have K, = h(uj — Q)L Z , 

L z = Em=o ma m M V' Since = °' the S round 

state of the Hamiltonian H — D,L Z can be determined 
once the spectrum of V as a function of L z (Vl) is ob- 
tained. 

Before proceeding, we define ^(r) = [— ^i(r) + 
l/Li(r)]/>/2, Ay(r) = ^_i(r)]/V2, A z (r) = 

^o(r)- Under a spin rotation, A rotations as a 3D vector 
in (xyz) space. Thus, states with total spin (5 = 0, 1, 
and 2) formed by two bosons at i"i and r 2 are : 

A\ x A, and A^A^ - \5i 3 A ' A' wnere A = A r i)- 

II. The role of internal degrees of freedom and 
the 'quenching' angular momentum L*: We have 
numerical diagonalized V up to N = 5 particles. The 
spectrum Vl of the N = 4 cluster with c 2 > is shown 
in figure 1 . It shows a clear correlation between spin and 
orbital angular momentum in the ground state. In the 
following, we shall explain the origin, the systematics, 
and the analytic structures of these states. 

Since c Q > |c 2 |, the minimum of the interaction Vij is 



zero, which occurs if the relative angular momentum I 
between bosons i and j is non-zero. For scalar bosons, 
Bose statistics demands £ > 2. A pair of spin-1 bosons, 
however, can make Vij — with t — 1 by making its spin 

part antisymmetric, i.e., forming a spin-1 state J(u\ — 
U2 ) A \ x A 1 0) ■ This shows that angular momentum states 
forbidden to scalar bosons will become accessible in the 
presence of internal degrees of freedom. Relative angular 
momentum, however, costs kinetic energy. One will then 
have to balance potential and kinetic energy. 

Among all the ground states, the state with minimum 
angular momentum L* that 'quenches' interaction (i.e. 
making V = 0) plays a special role. It will be referred 
to as the 'minimal quenching' state |$*). The vanishing 
of V can only occur if no two bosons can occupy the 
same point in space, meaning that B in cq. ([!]) contains a 
Laughlin factor W[z] = IlAr> i>:; > 1 (M i — Uj). For spinless 
particles, Bose statistics demand that W appears twice; 
resulting in a minimal quenching state B — W 2 , with 
L* = N(N — 1). For spin-1 bosons, the condition for 
quenching means 



B([u],[fi})=W[u]F([ u }M), 



(2) 



where F([u], [/i]) obeys Fermi statistics. Thus L* = 
N(N — l)/2 + P, where P is the lowest angular momen- 
tum possible for the fermionic wavefunction F. 

Since there are only three distinct spin states for spin- 
1 bosons, the only clusters whose quenching 'fermion' 
component F has L — are N = 2,3. For N — 2, we 
have |$*)( 2 ' = J u 12 A x A\®)i where My = ttj - Uj. 
It has (L* = 1,5 = 1). For N = 3, we have |$*)( 3 ) = 

Ju 12 u 23 u 31 A\ x A ■ A\°)> with ( L * = 3,5 = 0). For 
N > 3, the quenching fermionic function F must have 
non-zero angular momentum, for otherwise there will be 
two 'fermions' occupying the same state. One must there- 
fore populate particles in different spin and orbital states. 
This procedure gives rise to a unique minimal quenching 
structure, whose general form will be clear after we enu- 
merate a few cases. Writing = Jw D^lO) , we 
have 



(3) 



which gives (L* = 7,5 = 1), (L* = 12,5 = 1), (L* = 
18, 5 = 0), and (L* = 26, 5 = 1) for N = 4, 5, 6 and 7 
respectively. In contrast, the corresponding L* for scalar 
bosons are 12, 20, 30, 42 , about 100% higher. Since W 
is antisymmetric, only the antisymmetric part of D^ N ' 
will appear in The antisymmetrization turns 

the operator product into a determinant, just as fermion 
systems fl3|| . This prescription yields spin 5 = 0, 1, 1 and 
quenching angular momentum L* = 2A 2 /3 — N + S for 
A = 0,l,2(mod3). 
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Note that as particle number increases, more and more 
singlets made up of boson triplets appear in the minimal 
quenching state. This is a clear sign that the system 
becomes a spin liquid while evolving into a full fledged 
quantum Hall state. Note also that the Laughlin factor 
W forces all ground states with L > L* to have zero 
interaction energy. 

III. The route to quenching: We now turn to the 
ground states with L < L* . To illustrate the correla- 
tion between S and L, we denote the spin of the ground 
state of V with angular momentum L as Sl, and dis- 
play their value up to L* for a cluster of N bosons as 
[S L } N = (S ,Si,...,S L *) N . For c 2 > 0, we have (0, 1) 2 ; 

(1,1,1,0)3, (0,1,0,1, 0,1,1,1)4, (1,1,1, 1,1,1, 1,1,1, 

1,0,0, 1)5 ; and for c 2 < 0, we have (2, 1) 2 , (3,2,1,0) 3 , 
(4,3,2,1,2,1,1,1)4, (5,4,3, 3,1,1,3,1,1, 1,2,2, 1) 5 . 

To understand these structures, we note that since L < 
L* , there is not enough angular momentum to produce 
a Laughlin factor to prevent the coincidence of any two 
bosons. Still, it is desirable to generate as many boson 
pairs with relative angular momentum of unit strength as 
possible. This consideration then motivates some simple 
rules for possible candidates of the ground state, and 
for small clusters, they often narrow down to the exact 
structure. The rules are: (i) For given L, the ground state 
of interaction V will contain a maximum number of pairs 
with unit relative angular momentum, (ii) For c 2 > 0, (or 
c 2 < 0) the ground state of V will have a minimum (or 
maximum) total spin consistent with rule (i) and Bosc 
statistics. We have constructed states according to these 
rules for clusters up to N — 4 bosons. The results are 
given in Table 1, which reproduces the sequences [S'l]'- 4 -' 
listed above. When compared with the exact numerical 
results in the limit |c 2 |/c Q << 1, we found that most 
of states in Table 1 coincide with the exact numerical 
results, or otherwise have over 94% overlap with them. 
Since the interaction V is invariant under separate spin 
and orbital rotation, the correlation between S and L 
comes entirely from statistics, just as the Hund's first 
and second rule. 

In figure 2, we show the densities of a N — 5 clus- 
ters with L — 0, 3, 6, 9, 12, and the spin densities for the 
cluster with iV = 5,L = 3,S'=l. The emergence of 
a quantum Hall plateau in the density is clear near the 
quenching limit. (See also Section V). 

IV. The ground states as a function of 57: Ex- 
periments performed at fixed rotational frequency f2 are 
described by the Hamiltonian 7i = H — D,L Z . The ground 
state is given by the minimum of Hl = Vl + — fi)L, 
where Vl is the ground state energy of V with L z = L. 
If L were a continuous variable, and Vl a smooth curve, 
the optimum value L Q will satisfy h(u> — 17) = —OVl/QL, 
and d 2 VL/d 2 L > 0. In addition, we have dVt/dL — 
(<9 2 Vl/<9 2 L)l o . From figure 1, we see that aside from a 
small number of cusps (in this case at L — 5 and 7) the 
envelop of Vl is smooth and concave up, implying that as 
Q is increased the angular momentum increases in small 



(but discrete) steps, reaching L* at a critical frequency 
Q*. 

For Q > Q* the minimum of TLl within the interval 
(0 < L < L*) still occurs at L*. However, over the 
interval L > L* , we have Vl = 0, and hence TLl = ^(w — 
£l)L has a minimum at L*. This establishes Statement 
(D) in the opening, that for sufficiently rapid rotation 
the angular saturates at L*. 

V. The spin and quantum Hall signature of ro- 
tating clusters: Because of the large number (Q) of 
identical but phase incoherent clusters in the insulating 
phase, many quantities measured for the entire group of 
clusters coincide with the quantum mechanical averages 
of the same quantity within a single cluster, except that it 
is magnified by Q. This is best illustrated by considering 
the density of the entire group of clusters after the trap 
is turned off. Because of tight trapping, the expansion is 
mainly driven by the localization energy of the trap and is 
therefore ballistic. The field operator for the cluster (say 
at R = 0) then evolves as Vv( r ) = /dr'Q(r — r'; t)^(r'), 
where G(r;t) = U^i^at), G(x;t) = {M/2mMf/ 2 
exp[iMx 2 /2M]. Before the expansion, the field operator 

is ?/v(r) = J2e w e( r ^)fo( z ) a i^ anci tne single particle 
density matrix of the cluster is 

P»u(t) = (^(r)W(r)) = g^{r±)v{z), (4) 
9^(r±) = T,e \we(rx)\ 2 (a\ v a tlil ), (5) 

with v(z) — \fo(z)\ 2 , and rj^ = (x,y). In the in- 
sulating regime, the density matrix of the entire clus- 
ter collection is p^l(r) — J2RP^( r ~ where the 
sum is over all clusters. After the expansion, we have 
i>n(r,t) = J2i w d r ±-) fo( z ) a _Ui> w e(rx) = J G(x~x';t) 
G(y-y';t) w e (r' ± ) dr' ± , and J (z) = J 'G(z - z'; t) fo(z') 
dx'. This implies p^(r,t) ={ ipl(rt) ^(rt)) = g^Tx) 
v(z), where g M „(rjJ and v(z) are <? p „(rjJ and v(z) in 
eq.([3|) with wi and fo replaced by wi and f . Evaluat- 
ing g^y and v, we find that g^(r) = A(t) ~ 2 g^(r/A (t)), 

v(z) = A(t)- 1 i;(z/A(t)), where A(t) = y/l + {ujt) 2 . We 
then have 

p^(r,t) = A(t)- 3 E R P^([r-R]/A(t)) (6) 

The diagonal elements of eq.(|[) give the density of each 
spin component. Eq.([|) shows that the density after ex- 
pansion is simple a sum of all the expanded clusters, each 
of which obeys a scaling relationship. For times t >> 
1/w, each cluster has expanded to a size much larger than 
the dimension (A) of the origin cluster ensemble, eq.(^) is 
then well approximated by p^ l (r,t) = Q/(o;t) 3 p M (r/u>t) 
up to a correction of order O (A/auit). The density of 
each spin component of the entire expanded ensemble at 
long times therefore reproduces that of each individual 
cluster. (For magnetic field effects, see Q). 

The identity X^MUr) ( r /a) 2 i>^(r) = L z + N in the 
LLL implies that if L is the angular momentum of a 
cluster before expansion, then the mean square radius 
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at time t » 1/ui is r 2 = Jr 2 ± p all (r,t) = Qa 2 (ivt) 2 (L + 
N), and I 2 = Jz 2 j_p all (r,t) = Qa 2 (ujt) 2 N/2. One can 
therefore extract the angular momentum of the system 
from the shape of the entire expanded clusters. Moreover, 
the difference of mean squared radius between successive 
L states as increases is Qa^cot) 2 . 

We would like to thank Immanuel Bloch for discus- 
sions. This work is supported by NASA Grants NAG8- 
1441, NAG8-1765, and by NSF Grants DMR-0109255, 
DMR-0071630. 

Notes added: Concurrent with the submission of this 



paper, a preprint (cond-mat/0203061) by B. Paredes, P. 
Zoller, and I. Cirac has appeared. Apart from the gen- 
eral discussion of fermionization, there are little overlaps 
between our papers. 
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Table 1: The ground states constructed from the rules 
in Section III. The ground state is written as |<I>) = 

Jip L D^\0), using 0i2 =A X ■ A 2 , T 123 =A 1 kA 2 ■ A 3 , 



B 12 =A 1 x A 2 ,A+ = A X 



iA 



N>i>j> 



l(«i 



Uj), Uij = Ui — Uj. All states listed are exact eigenstates 
except those with L = 4, 5, 6. For L — 4, the last factor 
of ip can be either U24 or U34. The exact ground state 
is a combination of both. The state of L = 5 and 6 has 
99.6% and 94% overlap with the exact ground state. 
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If the expansion takes place in the presence of a magnetic 
field B, and if the original trap has cylindrical symmetric 



along z, [V = iM(u;ir : 
ized to pf{v,t) = 3-^3 



2 2 



eq.q 

U 



i|) will be general 



A 



\±At)= v / i + (^,^) 2 . u= e- l0 -"\ 6 = yBt, 

where 7 is the gyromagnetic ratio. At long times (ui±t >> 
1, 0J z t » 1), we have pf(r,t)= [Q / '(u) z w 2 j_t 3 )] L> MQ 
p a /3(r±/uj±t, z/uizt) Up„, which provides information 
about the density matrix of the pre-expansion cluster. 
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Figure 1: Vl for an N = 4 cluster with C2 > 0. 
Figure 2a: The column densities of a N = 5 cluster with 
with L = 0, 3, 6, 9, 12. The central density of these clus- 
ters decreases with increasing L. The L = 12 data is in 
bold. Figure 2b: The spin densities of the L — 3 cluster 
in fig.2a.; (S = 1). p = 1/na 2 . 
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